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Reynolds  Averaged  Navier  Stokes  simulations  have  been  performed  to  examine  modeling 
upgrades  for  scramjet  flowpath  predictions.  A  flush,  non-reacting  hydrogen  fuel  injector 
flowfield  was  used  as  the  model  problem,  and  an  LES  simulation  of  the  problem  was  used 
for  comparison  purposes.  Calculations  were  first  performed  examining  the  effect  of  Schmidt 
number  with  a  constant  Prandtl  number.  Next  the  effect  of  the  compressibility  correction 
used  was  examined.  These  findings  indicated  that  for  this  injector  configuration,  the  effect 
of  the  compressibility  had  a  major  impact  on  the  solution,  and  that  the  average  Schmidt 
number  of  about  0.45  compared  closely  to  the  LES  simulation  results.  Next,  a  new  scalar 
fluctuation  model  was  used  to  obtain  local  values  of  Prandtl  and  Schmidt  number  whose 
values  were  found  to  vary  significantly  across  the  fuel  jet  mixing  layer.  The  turbulent 
Prandtl  number  was  found  to  vary  between  0.4  to  0.9,  and  the  turbulent  Schmidt  number 
varied  from  0.6  to  1.2.  Finally,  a  comparison  was  performed  using  an  unstructured  flow 
solver  with  grid  adaptation.  This  technique  is  now  being  used  to  obtain  grid  resolved 
solutions  in  a  systematic  and  straightforward  manner  in  our  design  studies. 


Nomenclature 


E 

= 

Oxygen  Entrainment  Parameter 

K 

= 

Turbulent  Kinetic  Energy  (m2/s2) 

Kt 

= 

Temperature  Fluctuation 

Ky 

= 

Species  Fluctuation 

T 

= 

Temperature  (K) 

U,  V,  w 

= 

Velocity  (m/s) 

X 

= 

Streamwise  Coordinate  Direction 

Yj 

= 

Mass  Fraction  of  species  i 

5 

= 

Boundary  Layer  thickness  (0.1524  cm)  /  Grid  unit 

8,  St,  £y 

= 

Turbulent  Dissipation 

Prt 

= 

Turbulent  Prandtl  Number,  ratio  of  viscous  to  thermal  diffusion 

Sct 

= 

Turbulent  Schmidt  Number,  ratio  of  viscous  to  mass  diffusion 

Le 

= 

Turbulent  Lewis  Number,  ratio  of  mass  to  thermal  diffusion 

LES 

= 

Large  Eddy  Simulation 

RANS 

= 

Reynolds  Averaged  Navier  Stokes 
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TI.  Introduction 

his  paper  primarily  addresses  applying  scalar  fluctuation  modeling  to  a  non-reacting  hydrogen  fuel  injection 
problem  emulating  the  environment  of  a  hypersonic  scramjet  combustor.  A  companion  paper  1  discusses  details  of 
the  modeling  and  its  validation  for  some  unit  problems.  We  briefly  review  the  equations  implemented  in  Section  II. 
In  Section  III,  we  describe  a  detailed  LES  simulation  of  this  fuel  injection  problem  that  builds  upon  earlier  work  by 
Kannepalli  et  al.  2  and  uses  recycling/rescaling  methodology  to  interface  the  LES  jet  simulation  with  that  of  the 
approach  RANS  boundary  layer.  The  LES  simulation  is  computationally  intensive  and  has  progressed  to  a  state 
where  RMS  mean  flow  quantities  are  converged  but  turbulent  stresses  and  scalar  variances  are  still  somewhat  noisy, 
but  are  yielding  reasonable  values.  Extracting  local  values  of  Prt,  Sct  and  Le  from  the  LES  simulation  entails 
invoking  a  gradient  transport  hypotheses,  i.e.,  aa  (a  =  Prt  or  Sct  where  a  is  the  temperature  or  species  mass 
fraction  respectively)  is  obtained  from 


u'a'  dui  / dxi 
u'u'j  da/dxi 


(1) 


where  division  of  second-order  correlations  is  entailed.  Thus,  aa  can  be  very  noisy  and  the  LES  calculation  must  be 
run  further  to  obtain  “smooth”  values.  This  is  now  in  progress. 

RANS  simulations  were  performed  on  this  same  fuel  injection  geometry  with  a  similar  boundary  layer  inflow 
profile.  Calculations  were  performed  first  looking  at  the  effect  of  Lewis  number,  using  constant  values  of  Prt  and 
Sct,  and  then,  using  the  scalar  fluctuation  equations  to  provide  local  values.  Calculations  were  performed  with  and 
without  a  compressibility-correction  [3]  activated  comparing  the  RANS  calculations  with  the  averaged  LES  data. 
These  comparisons  are  discussed  in  Section  IV.  Our  findings  are  as  yet  inconclusive,  but  suggest  that  modifications 
to  the  basic  compressibility  correction  are  needed  for  this  class  of  flows. 

Finally,  in  Section  V,  we  show  applicability  of  our  unstructured  code,  CRUNCH®,  to  this  problem.  For  realistic 
applications,  CRUNCH®  has  been  used  in  the  fuel  injection  zones  of  various  scramjet  concepts,  since  the  geometry 
of  multiple  fuel  injection  ports  is  more  readily  modeled  using  an  unstructured  approach.  Using  the  unstructured 
solver  also  allows  the  advantage  of  our  grid  adaptation  capabilities,  which  have  been  found  to  have  a  first  order 
effect  on  parameters  such  as  combustion  efficiency,  and  thus  thrust 3. 

II.  Problem  Definition  and  Computational  Methodology 

In  this  section  a  description  of  the  numerical  problem  will  be  discussed.  The  first  part  will  describe  the  fuel 
injection  geometry  that  was  used  to  obtain  the  numerical  results.  The  part  section  will  describe  the  initial  and 
boundary  conditions  for  the  domain.  The  last  part  will  describe  the  methodology  of  the  computational  models  that 
were  used. 

A.  Geometry 

Figure  1  shows  the  fuel  injection  region  and  channel  grid  structure.  The  grid  is  scaled  based  on  5,  the  thickness 
of  the  incoming  turbulent  boundary  layer.  The  channel  is  358  long,  108  high,  and  58  wide.  A  flush  wall  injector, 
angled  at  25  degrees  relative  to  the  channel  wall,  is  placed  58  downstream  of  the  inflow  plane  along  the  channel 
centerline.  The  injector  is  a  square  of  0.58  on  each  side.  The  symmetry  plane  down  the  centerline  of  the  injector  was 
utilized  so  that  only  half  of  the  domain  needed  to  be  modeled.  The  grid  was  constructed  as  one  block  of  2.8  million 
cells,  with  I-blanking.  The  grid  was  clustered  near  the  wall  and  where  the  injector  flowed  into  the  channel. 

This  same  general  geometry  and  gridding  was  used  for  both  the  RANS  and  LES  calculations.  The  LES 
calculation  does  not  use  the  symmetry  plane  down  the  middle  of  the  channel,  and  needed  to  model  the  entire  channel 
width  of  108.  Also,  due  to  the  expense  of  the  LES  calculation  compared  to  the  RANS,  only  158  downstream  of  the 
injector  were  modeled.  The  LES  calculation  also  included  a  158  long  and  28  wide  extension  upstream  on  the 
channel  that  was  used  to  generate  the  unsteady  turbulent  boundary  layer. 
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Figure  1.  Numerical  Model  Geometry 


B.  Inflow  and  Boundary  Conditions 

The  inflow  and  injector  conditions  are  given  in  Table  1.  The  channel  inflow  was  initialized  as  air  composed  of 
02  and  N2,  with  the  species  mass  fractions  of  0.2323  and  0.7677  respectively.  A  turbulent  boundary  layer  of 
thickness  one  8  was  created  using  a  l/7th  power  law  distribution  for  the  velocity,  and  a  temperature  profile  using  the 
Crocco-Busemann  relationship  5.  The  incoming  Mach  number,  temperature,  and  pressure  are  representative  of  the 
conditions  for  the  inflow  to  a  scramjet  combustor.  The  injector  conditions  were  set  for  a  pure  hydrogen  (H2)  inflow 
at  near  sonic  conditions.  The  channel  wall  was  a  no- slip,  isothermal  surface  having  a  prescribed  temperature  of 
700K.  The  y+  values  ranged  from  0.2  to  0.4  along  the  wall  surface.  The  channel  exit  was  a  standard  extrapolated 
supersonic  outflow. 


Table  1.  Inflow  and  Injector  Conditions 


Inflow 

Injector 

u  (m/s) 

3131.7 

1920.0 

Mach 

4.5 

1.05 

Ts  (K) 

1200.0 

575.0 

Ps  (psi) 

14.7 

325.0 

Y02 

0.2323 

0 

yN2 

0.7677 

0 

Yh2 

0 

1 

C.  Computational  Methodology 

For  this  study  two  CFD  codes  were  used,  one  using  a  structured  methodology,  and  the  other  using  an 
unstructured  methodology.  The  structured  code  was  the  CRAFT  CFD®  solver  4,  with  Table  2  listing  the  codes 
features.  This  code  was  used  for  both  the  RANS  and  LES  calculations.  For  the  RANS  simulations  a  k-c  turbulence 
model,  shown  in  Equation  2  and  3,  was  used  with  the  So-Zhang-Speziale  6  near  wall  model.  A  compressibility 
correction  was  also  available  based  on  a  modification  to  the  Sarkar  model  1 .  The  flow  was  modeled  as  a  non¬ 
reacting  mixture  of  thermally  perfect  gases  (02,  N2,  and  H2). 


3 

American  Institute  of  Aeronautics  and  Astronautics 


AIAA-2005 -0432 


Table  2.  Relevant  Features  of  CRAFT  CFD®  Navier-Stokes  Code 


NUMERICS/ 

PARALLEL 

PROCESSING 

•  1D/2D/AXI/3D  Finite- Volume  Discretization 

•  Implicit,  ADI  and  L/U,  Upwind  (Roe/TVD)  Formulations 

•  Fully  Implicit  Source  Terms/Boundary  Conditions 

•  PNS  Run  Option 

•  Domain-Decomposition  Parallel  Architecture  with  MPI 

•  Preconditioning  Extensions 

GRID  FEATURES 

•  Grid  Patching/Blanking  for  Complex  Geometries 

•  Noncontiguous  Grid  Interfacing  with  Flux  Preservation  Across  Domains 

THERMO¬ 

CHEMISTRY 

•  Real  Gas  Mixtures  (Calorically  and  Thermally  Imperfect/ J  ANN  AF  Thermo  Tables) 

•  Generalized  Finite-Rate  Chemistry  and  One-step  Option 

•  Fully  Implicit  Source  Term  Linearization 

•  Probability  Density  Function  Turbulent  Chemistry  Modeling 

•  Stiff  Chemistry  Solver 

TURBULENCE 

•  k-c  /  EASM  Nearwall  Functionality 

•  Unified  Compressibility  and  Vortex  Stretching  Corrections 

•  Simplified  and  Advanced  Transition  Modeling 

•  Particle  Dispersion  Formulations 

•  Hybrid  RANS/LES 

MULTIPHASE 

FLOW 

•  Nonequilibrium  Particle/Droplet  Solvers  (Eulerian  and  Lagrangian  Formulations) 

HIGH  ALTITUDE 
MODELING 

•  Vibrational  Non-equilibrium  Modeling 

•  Slip  Wall  Boundary  Conditions 

dpk_  dPujk 

dt  dx . 

J 

d£_  a?uje 

dt  dx . 

J 


d 

V 

p  H - 

dk 

dx  . 

dx  . 

(7  j 

J 

{  k) 

J  \ 

d 

V 

d£ 

dx  . 

p  H - 

G 

dx  . 

J 

\  £ ) 

J 

+pk 


+p 

£ 


-D 

£ 


+  SSj 

k 


+  55 
£ 


(2) 


(3) 


To  improve  the  capability  for  predicting  scalar  transport  parameters  in  complex  geometries,  an  earlier 
temperature  variance  predictive  methodology  8  which  operated  within  the  framework  of  a  RANS-based  k-c 
turbulence  model,  was  extended  by  Brinkman  The  extended  model  solves  for  both  the  temperature  and  species 
variance,  from  which  the  scalar  transport  parameters  Prt  and  Sct  are  derived.  Transport  equations  for  the  temperature 
variance  kT  and  its  dissipation  rate  cT  were  solved,  coupled  to  the  k-c  turbulence  model.  The  transport  equations  for 
kT  and  sT  are: 


d(pkT)  d(pujkT ) 

dt  dx . 

J 


d 


dx  . 
J 


i  \ 

^  °t  i 

dkT 

OL  -| - 

(Jk,T  J 

dx  . 

J  . 

+  K,  ~Dk 


(KP£t} 

dt 


d(p  Uj£t  ) 


dx  . 
J 


dx  . 
J 


at 

d£T 

p 

ceH - - 

dx  . 

J 

O  ej  j 

RR 


D-r  +  Ck+^eT 


(4) 


(5) 
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where  the  compressibility  correction  accounted  for  in  the  k-s  model  is  carried  over  with  the  expression  Ck  and  ^sT  is 
a  near-wall  damping  function  included  to  capture  low-Reynolds  number  behavior.  The  turbulent  Prandtl  number 
with  near- wall  damping  is  calculated  as 


f n  k  Sj 

c\f\iek 


where  fM  and  fa  are  the  near- wall  damping  term  for  turbulent  momentum  diffusivity  and  temperature  diffusivity 
respectively,  and  CA  and  are  proportionality  constants.  A  similar  formulation  is  implemented  for  the  species 
mass-fraction  variance  kY  and  its  dissipation  rate  sY ,  to  calculate  a  turbulent  Schmidt  number,  Sct,  and  mass 
diffusivity.  For  a  detailed  explanation  and  definitions  of  all  the  variables  in  equations  5  and  6,  see  Reference  [1]. 

For  LES  applications,  the  CRAFT  CFD®  1  solver  was  implemented  with  an  upwind-biased,  Roe-flux- 
extrapolation  procedure  that  has  been  extended  to  fifth  order 9,10  for  the  inviscid  and  fourth  order  central  differencing 
for  the  viscous  terms.  Temporally,  the  code  includes  both  a  fourth  order  Runga-Kutta  scheme  and  a  second-order 
three-factor  Approximate  Factorization  (AF)  implicit  scheme.  For  the  studies  described,  the  AF  scheme  was  used 
and  applied  with  sub-iterations  to  remove  the  splitting  error.  For  subgrid  modeling,  the  code  includes  a 
compressible  version  of  the  algebraic  Smagorinsky  model  11  as  well  as  a  one-equation  model  of  Menon  12 .  The  one 
equation  model  solves  a  transport  equation  for  subgrid  turbulent  kinetic  energy,  ksgs .  The  subgrid-scale  stresses  are 
then  modeled  using  an  eddy  viscosity  approach  based  on  k?8S.  In  order  to  stabilize  the  code  in  the  vicinity  of  strong 
gradients,  such  as  shock  waves,  a  modification  of  the  classic  Jameson,  et  al.13  2-4  dissipation  scheme  was  used  with 
a  cutoff  switch  that  is  tuned  to  add  dissipation  only  at  shocks  and  discontinuities2. 

For  the  unstructured  simulations,  the  CRUNCH  CFD®  4  solver  was  used.  Table  3  lists  these  code  features. 
CRUNCH  CFD®  is  a  hybrid  unstructured  RANS  solver  employing  a  k-c  turbulence  model  similar  to  the  CRAFT 
CFD®  solver.  The  compressibility  correction  is  based  on  an  unmodified  Sarkar  model7. 


Table  3.  Current  Features  of  the  CRUNCH  CFD®  Code 


NUMERICS 

•  Finite- Volume  Roe/TVD  Flux  Construction,  Vertex  Storage 

INTEGRATION 

•  Explicit  Four-Step  Runge-Kutta,  Implicit  GMRES,  Gauss-Seidel 

GRID  ELEMENTS 

•  Tetrahedral,  Hexahedral,  Prismatic,  Pyramid 

PARALLEL  PROCESSING 
CAPABILITIES 

•  Domain  Decomposition  MPI,  Independent  Grids  with  Noncontiguous 
Interfacing,  Automated  Load  Balancing 

DYNAMIC  GRID 
CAPABILITIES 

•  Node  Movement  Solver  (Implicit  Elasticity  Approach),  Automated 

Embedding,  Sliding  Interfaces 

GRID  ADAPTION 

•  Variable  Element  Grid  Refinement  using  Delaunay  and  cell  subdivision 

Procedures,  Automated  Load  Balancing  of  Adapted  Grid 

THERMOCHEMISTRY 

•  Multi-component  Real  Gas  Mixtures,  Finite-Rate  Kinetics 

•  kc  /EASM  Formulations  with  Compressibility/Vortical  Upgrades 

TURBULENCE  RANS/LES 

•  LES  Subgrid  Scale  Models  -  Algebraic  and  One-equation 

•  Unified  Hybrid  LES/RANS  Framework 

MULTIPHASE  FLOW 

•  Nonequilibrium  Particle/Droplet  Solvers  (Eulerian  and  Lagrangian 
Formulations),  Cavitation  Model  with  Bubble  Dynamics 

III.  LES  Simulation 

The  LES  version  of  the  CRAFT  CFD®  solver  was  used  to  solve  the  3D  fuel  jet  injection  problem.  To  obtain  a 
turbulent  and  unsteady  boundary  layer  inflow  for  the  calculation,  a  computational  domain  upstream  of  the  main 
channel  was  used.  A  domain  length  of  155  in  the  streamwise  direction  was  required  to  set  up  the  unsteady 
supersonic  boundary  layer  via  the  recycling/rescaling  technique  2.  A  wall-normal  dimension  of  105  and  spanwise 
dimension  of  25  (with  periodic  boundary  conditions)  was  sufficient  to  generate  a  “bona-fide”  turbulent  boundary 
layer. 
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Figure  2  shows  a  side  view  along  the  length  of  the  channel,  taken  down  the  centerline  of  the  jet.  This  figure 
shows  the  representative  distances  at  which  axial  cuts  of  the  domain  are  taken  for  comparison  purposes.  The  x 
coordinate  locations  shown  in  the  figure  are  relative  to  the  start  of  the  jet  and  are  in  grid  units,  where  1  grid  unit,  8, 
is  representative  of  the  thickness  of  the  incoming  boundary  layer.  The  coordinate  system  origin,  x  =  0,  begins  where 
the  injector  initially  makes  contact  with  the  channel  wall.  Since  the  injector  is  at  a  25  degree  angle  with  respect  to 
the  channel  wall,  the  end  of  the  injector  is  at  about  x  =  1.5.  The  channel  locations  selected  for  presentation  are 
located  at  18,  38,  58,  and  78.  Since  the  injector  was  designed  to  have  a  diameter  of  14  8,  these  distances  also 
represent  2,  6,  10,  and  14  jet  diameters  downstream. 


Figure  3  shows  the  contours  for  the  mean  values  of  velocity,  u/uoo,  temperature,  T/T^,  hydrogen  mass  fraction, 
YH2,  turbulent  kinetic  energy,  Vk/  Uoo,  temperature  fluctuation,  Vkt/  T^,  and  species  fluctuation,  Vky,  at  the  cutting 
planes  defined  in  figure  2.  The  turbulent  kinetic  energy,  temperature  and  species  fluctuations  are  calculated  from  the 
LES  simulation  as: 


K  =  y2(u'u'  +  v'v'  +  w'w') 

(7) 

Kt=T'T' 

(8) 

II 

(9) 

Since  the  LES  calculation  modeled  the  entire  domain,  asymmetry  is  seen  around  the  jet  centerline.  These  figures 
show  the  characteristic  vortex  formed  by  a  jet  entering  a  supersonic  flow.  This  vortex  assists  in  the  mixing  of  the 
fuel  jet  with  the  airflow  through  the  channel.  An  additional  component  of  the  mixing  is  due  to  the  mass  diffusion  of 
the  fuel  into  the  air.  Figure  4  shows  a  three-dimensional  view  of  the  LES  simulation.  The  iso-surfaces  of 
temperature  (1500  K)  are  in  green,  density  (0.5  kg/m3)  is  in  orange  and  the  mass  fraction  of  hydrogen  (0.5)  is  in 
white.  The  density  iso- surfaces  show  the  spherical  shock  wave  structure  propagating  downstream  of  the  jet.  This 
shock  wave  fully  encloses  the  injector  jet.  This  figure  shows  the  asymmetry  associated  with  LES  simulations. 

IV.  RANS  Simulation  with  the  CRAFT  CFD®  Solver 

The  results  presented  in  this  section  are  for  the  structured  RANS  solver.  Two  subsections  will  be  discussed.  The 
first  set  of  calculations  look  at  examining  the  effect  of  varying  the  Schmidt  number  while  maintaining  a  constant 
Prandtl  number.  Also  looked  at  in  this  subsection  is  the  effect  of  the  compressibility  correction  on  the  solution.  The 
second  subsection  will  examine  the  use  of  the  variable  Prandtl  and  Schmidt  number  model  discussed  above,  with 
and  without  the  compressibility  correction.  Both  of  these  sets  of  calculations  will  be  compared  back  to  the  LES 
simulation,  which  will  be  used  like  a  baseline  calculation. 
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(e).  VKt/T^  (f).  Vky 


Figure  3.  LES  Profile  Contours. 
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Figure  4.  LES  Iso-Surfaces. 


A.  Constant  Prandtl  Number  Calculations 

The  calculations  listed  in  Table  4  performed  with  specified  values  of  transport  parameters  for  this  study  as  Cases 
1-4.  The  turbulent  Prandtl  was  chosen  for  all  cases  to  be  0.9,  a  value  that  is  commonly  used  for  boundary  layer 
applications.  The  Schmidt  number  varied  between  0.9  and  0.45,  making  the  Lewis  number  1.0  or  2.0  respectively. 
Lowering  the  turbulent  Schmidt  number  (raising  the  Lewis  number  since  Le  =  Prt/Sct )  permits  the  species  to  diffuse 
faster  than  the  temperature,  which  for  burning  cases  can  result  in  earlier  burning,  see  for  example  Baurle  and 
Eklund14.  The  Schmidt  number  (or  Lewis  number)  is  the  parameter  controlling  turbulent  diffusion  of  the  fuel  jet  into 
the  air  stream.  The  compressibility  correction  was  also  a  factor  considered  in  this  study.  This  correction  suppresses 
the  turbulent  mixing  rate  and  reduces  the  entrainment  of  air  into  the  jet/air  mixing  layer. 


Table  4.  Inflow  and  Injector  Conditions 


Case 

Compressibility 

Pr, 

Sct 

Le 

1 

ON 

0.90 

0.90 

1.0 

2 

ON 

0.90 

0.45 

2.0 

3 

OFF 

0.90 

0.90 

1.0 

4 

OFF 

0.90 

0.45 

2.0 

5 

ON 

Variable 

Variable 

Variable 

6 

OFF 

Variable 

Variable 

Variable 

Figure  5  shows  a  comparison  of  the  oxygen  entrainment  shown  in  into  the  fuel  jet  mixing  layer  for  these  four 
constant  Prandtl  number  cases  as  compared  to  the  value  from  the  LES  simulation.  The  entrainment  parameter,  E,  is 
defined  in  Eqn.  10  as  the  mass  flux  average  of  the  amount  of  02  that  is  mixed  into  the  fuel  jet.  The  fuel  jet  is  defined 
as  the  computational  domain  where  the  mass  fraction  of  fuel  is  greater  than  0.5%. 


X 

Figure  5.  Entrainment  of  02  Into  the  Jet. 
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Comparing  Cases  1  and  2,  which  have  compressibility  on,  with  the  Lewis  number  varying  from  1  to  2,  the 
entrainment  into  the  jet  is  shown  to  increase  due  to  the  increase  in  mass  diffusion,  which  is  keyed  to  the  value  of  the 
eddy  viscosity  divided  by  the  Schmidt  number.  Comparing  Cases  1  to  3,  which  have  the  same  Lewis  number,  but 
the  compressibility  correction  is  not  used  for  Case  3,  the  entrainment  increases  because  the  overall  mixing  is  faster 
since  the  eddy  viscosity  itself  is  larger.  Examining  the  entrainment  between  Cases  2  and  3  indicates  that  both  faster 
mixing  and  increased  species  diffusivity  increase  the  entrainment  by  the  same  magnitude,  since  the  values  of  the 
turbulent  diffusivity  is  the  same.  However,  in  reacting  cases,  the  effects  on  the  burning  may  not  be  the  same.  Based 
on  the  entrainment  parameter,  the  Case  4,  which  has  a  Schmidt  number  of  0.45  and  the  compressibility  correction 
off,  is  seen  to  provide  mixing  that  is  comparable  to  the  LES  simulation.  For  this  injector  configuration,  the  LES 
simulation  is  indicating  that  compressibility  has  a  major  effect,  and  that  the  lower  Schmidt  number  produces  results 
more  in  accord  with  the  LES  simulation. 


Figure  6  shows  a  comparison  of  the  four  cases  at  the  x=5  cutting  plane  for  the  velocity,  temperature,  mass 
fraction  of  fuel,  and  turbulent  kinetic  energy.  This  figure  should  be  compared  to  Figure  3  a-d  for  the  LES  simulation. 
Case#l  Case  #2  Case#l  Case  #2 


Case  #1 


Case  #2 


Case  #1 


Case  #2 


(c).  Yh2  (d).  Vk/Uq0 

Figure  6.  Contours  at  X=5  for  the  Constant  Prandtl  Number  Simulations. 
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These  figures  all  resemble  the  LES  contours,  but  the  turbulent  kinetic  energy  levels  for  the  cases  with  the 
compressibility  correction  off  compare  better  to  the  LES  calculation  than  the  cases  with  the  compressibility 
correction  on.  The  contours  of  H2  mass  fraction,  Ym,  also  show  that  the  Case  4  closely  matches  the  LES  cutting 
plane  contours.  Comparing  the  temperature  profiles  for  all  the  cases  to  the  LES,  it  is  seen  that  the  LES  has  higher 
temperature  values  in  the  boundary  layer  and  shock  structure  surrounding  the  jet.  This  was  caused  by  the  slightly 
different  inflow  profiles  being  used  for  the  RANS  and  LES  simulations.  The  RANS  used  the  l/7th  power  law  for 
velocity,  and  Crocco-Busemann  for  the  temperature  profile,  while  the  LES  used  a  recycling/rescaling  technique. 


B.  Variable  Prandtl  /  Schmidt  Number  Model 

As  described  above,  the  variable  Prandtl  number  model  solves  for  the  temperature  variance  and  dissipation.  The 
variable  Schmidt  number  model  additionally  solves  for  the  species  variance  and  dissipation.  From  these  two  sets  of 
quantities,  the  Prandtl  and  Schmidt  numbers  can  be  calculated.  Figure  7  shows  a  comparison  between  the  2D 
simulation  of  this  case  by  Brinkman  1  (transverse  slot  jet,  rather  than  a  discrete  jet  from  a  square  duct)  and  the  3D 
simulation  presented  here  for  the  normalized 
temperature  fluctuation.  For  the  2D  case,  all  the 
mixing  occurs  in  the  shear  layer  between  the  hot  air 
stream  and  the  cooler  fuel,  and  in  the  separation  zone 
that  forms  upstream  of  the  jet.  These  areas  show  high 
values  of  the  temperature  fluctuation.  Between  the 
shear  layer  and  the  channel  wall,  there  is  no 
mechanism  in  the  2D  flow  to  mix  the  fuel  and  air,  so 
little  to  no  temperature  fluctuation  is  seen.  The  3D 
simulation,  taken  down  the  centerline  of  the  jet,  shows 
the  high  values  as  the  jet  initially  flows  in  the  channel. 

The  horseshoe  vortex  provides  a  mixing  mechanism 
for  the  jet,  which  allows  fluctuations  across  the  fuel 
jet.  Allowing  for  the  differences  in  the  2D  versus  3D 
nature  of  the  flowfield,  these  simulations  compare 
well  to  each  other  in  the  general  structure  of  the 
temperature  fluctuation  field.  However,  the  3D  has 
greatly  enhanced  turbulent  kinetic  energy  due  to  the 
vortical  mixing  which  alters  all  the  fluctuation  levels  and  flow  structure. 

Table  4  lists  the  two  Cases,  5  and  6,  which  were  performed  for  this  study.  Both  of  these  cases  use  the  variable 
Prandtl  and  Schmidt  number  models,  but  Case  5  has  the  compressibility  correction  on  and  Case  6  does  not.  Figure  8 
shows  the  oxygen  entrainment  for  these  cases,  compared  to  the  LES  and  to  Case  1 ,  which  had  compressibility,  and  a 
constant  Lewis  number  of  1.  The  trend  seen  in  the  previous  entrainment  plot,  that  turning  the  compressibility 
correction  off  increases  the  mixing,  is  seen  in  this  plot  as  well.  The  magnitude  of  the  increase  is  not  as  great,  and 
could  be  due  to  the  coupling  of  the  turbulence  through  the  scalar  variance  equations. 

0.23 


Figure  7.  Temperature  Fluctuation  Comparison 
between  2D  (top)  and  3D  (bottom)  Models. 


0.17 


0.11 


Figure  8.  Entrainment  for  Variable  Prandtl  Number  Model. 
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Figure  9  shows  a  comparison  of  the  x  =  5  cutting  plane  between  the  LES  simulation  and  Case  6.  Shown  are  the 
velocity  and  turbulent  kinetic  energy,  temperature  and  temperature  fluctuation,  and  fuel  species  and  species 
fluctuation.  Case  6  was  chosen  over  Case  5  since  it  was  shown  in  the  previous  section  that  the  cases  with  the 
compressibility  correction  not  active  compare  better  in  turbulent  kinetic  energy  to  the  LES  simulation.  This  can  be 
seen  when  examining  Figure  9a.  Figure  9b  shows  the  temperature  and  temperature  fluctuation.  Both  of  these 
quantities  compare  well  to  the  LES  simulation,  with  the  LES  simulation  having  slightly  higher  values  for  the 
fluctuation.  The  fuel  species  and  fluctuation  in  Figure  9c  also  compare  well  in  shape.  The  LES  simulation  has  about 
double  the  value  of  the  fluctuation  compared  to  the  RANS  calculation.  Since  the  Schmidt  number  is  proportional  to 
the  square  root  of  the  reciprocal  of  the  species  fluctuation,  the  higher  species  fluctuation  in  the  LES  simulation 
would  translate  into  a  lower  Schmidt  number.  The  lower  Schmidt  number  causes  more  species  diffusion,  which  is 
represented  in  figure  9c  by  the  less  red  area  of  the  species  mass  fraction. 


(c) 

Figure  9.  Comparison  of  the  RANS  Variable  Prandtl  /  Schmidt  Number  Model  to  LES. 

Figures  10  and  1 1  show  the  jet  centerline  view  and  x  =  5  cutting  plane  view  of  the  Prandtl,  Schmidt,  and  Lewis 
numbers.  The  Prandtl  number  varies  between  0.4  -  0.9  in  the  fuel  jet  region.  The  Schmidt  number  also  varied  in  the 
jet  region,  with  values  ranging  from  1.2  -  0.6.  The  lower  values  occurred  in  the  top  shear  region  of  the  jet,  with 
higher  values  around  the  interior  of  the  jet.  The  core  region  of  the  vortex  had  low  values.  The  regions  where  the 
most  mixing  occurred  had  the  lowest  values,  since  this  is  where  the  species  fluctuation  would  be  greatest.  The  Lewis 
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number,  which  is  Prt  /  Sct  has  values  near 
2  in  the  outer  regions  of  the  fuel  jet  close 
to  the  injection  point,  and  these  high 
values  decrease  as  the  jet  propagates 
downstream.  The  interior  of  the  jet  has 
values  that  range  from  0.5  to  1.0,  with  a 
large  portion  of  the  jet  being  in  the  1.0 
region. 

These  studies  are  preliminary.  The 
scalar  fluctuation  model,  calibrated  by 
very  basic  unit  problem  data  sets1,  is 
certainly  operating  in  a  qualitatively 
proper  manner  for  a  3D  fuel  jet  problem 
of  this  complexity.  The  LES  calculation 
itself  can  be  taken  as  “correct”,  but  is 
still  not  entirely  converged. 


0.00  0.50  1.00  1.50  2.00 


Figure  10.  Variable  Prandtl  /  Schmidt  Number  Model  Centerline 

Profiles. 
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V.  RANS  Simulations  using  CRUNCH  CFD®  Solver 

The  calculations  performed  in  the  above  section  were  all  done  using  a  structured  grid,  RANS  solver,  CRAFT 
CFD®.  In  this  section,  results  obtained  using  an  unstructured  solver  CRUNCH  CFD®  will  be  exhibited.  In  many 
scramjet  and  fuel  injector  problems  in  general,  the  use  of  unstructured  grids  reduces  the  time  and  effort  required  to 
both  grid  certain  geometries,  and,  more  importantly,  to  obtain  grid  resolved  solutions  in  a  systematic  and 
straightforward  manner.  For  this  problem,  the  use  of  a  structured  grid  with  one  grid  block  required  I-blanked 
regions.  The  same  grid  structure  was  used  in  the  unstructured  solver,  but  the  I-blanked  regions  were  not  required, 
which  reduced  the  grid  size.  The  boundary  layer  inflow  and  jet  inlet  conditions  between  the  two  codes  were  nearly 
identical.  Both  codes  had  the  Prandtl  number  constant  at  0.9,  a  constant  Lewis  number  at  1.0,  and  the 
compressibility  correction  activated.  Figure  12  shows  a  comparison  between  the  CRAFT  CFD®  and  CRUNCH 
CFD®  solvers.  The  vortex  structure  is  slightly  more  bent  about  the  centerline  for  the  CRUNCH  CFD®  simulation, 
which  can  be  seen  in  the  velocity  and  species  mass  fraction  contour  plots.  Overall,  the  two  codes  compared  well  to 
each  other. 


Unstructured  Structured 


Figure  12.  Comparison  of  the  Unstructured  (left)  Solver  to  the  Structured  Solver  (right)  x  =  5.0. 


VI.  Conclusions  and  Future  Direction 

In  this  study  RANS  simulations  were  performed  and  compared  to  an  LES  simulation  on  the  same  geometry.  The 
LES  simulation  was  computationally  intensive  and  has  progressed  to  a  state  where  RMS  mean  flow  quantities  are 
converged  but  turbulent  stresses  and  scalar  variances  are  still  somewhat  noisy,  but  were  yielding  reasonable  values. 
The  first  set  of  calculations  performed  examined  the  effects  of  varying  the  Lewis  number  while  keeping  the  Prandtl 
number  constant.  These  results  showed  that  increasing  the  Lewis  number  resulted  in  greater  mixing,  but  not  to  the 
levels  of  the  LES  simulation. 

The  next  set  of  calculations  performed  examined  the  compressibility  correction.  It  was  determined  that  not 
having  the  correction  active  obtained  results  that  were  better  than  having  the  correction  active.  According  to  Spaid 
and  Zukoski  15  for  air  injection  into  an  air  stream,  the  compressibility  correction  is  required.  A  variation  of  this 
problem  with  Helium  injection  into  an  air  stream  1  showed  that  the  compressibility  correction  may  not  be  adequate 
to  capture  the  effects  of  the  large  air/Helium  density  gradient.  One  conclusion  from  this  is  that  the  compressibility 
correction  may  require  a  density  modification  to  account  for  the  interaction  of  different  gases. 

After  these  calculations  were  done,  a  new  variable  Prandtl  and  variable  Schmidt  number  model,  developed  by 
was  Brinkman  1  was  used.  This  model  showed  that  the  Prandtl  number  varied  from  0.4  to  0.9  in  the  jet  region.  The 
Schmidt  number  varied  from  0.6  to  1.2,  which  created  a  Lewis  number  variation  from  0.5  to  2.0  across  the  jet.  The 
modeled  temperature  variance  compared  well  in  general  shape  and  magnitude  to  the  LES  simulation.  The  species 
variance  compared  well  in  general  shape,  but  had  values  that  were  about  half  those  from  the  LES. 

The  final  comparison  was  performed  between  the  structured  solver  and  unstructured  solver  on  the  same 
computational  mesh.  Both  codes  were  run  with  the  same  Prandtl  and  Lewis  numbers,  but  had  slightly  different 
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version  of  the  compressibility  correction.  The  comparison  of  mean  flow  quantities  between  the  solvers  is  good. 
There  is  a  slight  difference  in  the  shapes  of  the  jets,  with  the  unstructured  solver  having  a  more  curved  jet. 
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